Generating Bounds for the Ground State Energy of the 
Infinite Quantum Lens Potential 



Carlos R. Handy 

Department of Physics and 
Center for Theoretical Studies of Physical Systems 
Clark Atlanta University 
Atlanta, Georgia 30314 

C. Trallero-Giner and Arezky H. Rodriguez 

Department of Theoretical Physics 
University of Havana 
10400 Havana, Cuba 



1 



Abstract 

Moment based methods have produced efficient multiscale quantization 
algorithms for solving singular perturbation/strong coupling problems. One 
of these, the Eigenvalue Moment Method (EMM), developed by Handy et 
al (Phys. Rev. Lett.55, 931 (1985); ibid, 60, 253 (1988b)), generates con- 
verging lower and upper bounds to a specific discrete state energy, once the 
signature property of the associated wavefunction is known. This method is 
particularly effective for multidimensional, bosonic ground state problems, 
since the corresponding wavefunction must be of uniform signature, and 
can be taken to be positive. Despite this, the vast majority of problems 
studied have been on unbounded domains. The important problem of an 
electron in an infinite quantum lens potential defines a challenging extension 
of EMM to systems defined on a compact domain. We investigate this here, 
and introduce novel modifications to the conventional EMM formalism that 
facilitate its adaptability to the required boundary conditions. 
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I. Introduction 



Self-assembled quantum dots (QDs), obtained by interrupted growth in 
strained semiconductors, offer an attractive and fascinating array of physi- 
cal properties (Leonard et al (1993,1994)). Differential capacitance (Drexler 
et al (1994), Miller et al (1997)), magnetic-conductance (Medeiros-Ribeiro 
et al (1997)), and optical experiments (Fafard et al (1994), Lee et al (2000)) 
demonstrate that electronic states are strongly confined inside such struc- 
tures. 

Typically, a lens geometry is assumed (Leonard et al (1994)), with 
a circular cross section of maximum radius a, and maximum thickness b; 
wherein, the charge carriers are confined by a hard wall (infinite) potential. 
The mathematical characterization of the energy levels of such nanostruc- 
tures is a delicate problem, particularly in the thin lens limit ^ — ^ 0, which 
corresponds to a singular perturbation regime. 

Recently, conformal analysis methods were used to solve the infinite 
quantum lens potential (Rodriguez et al (2001)). Preliminary results un- 
derscore the delicate nature of the thin lens regime. In order to better assess 
the accuracy of such methods, we have developed an eigenenergy bound- 
ing procedure that, at low order, yields exceptionally tight bounds to the 
discrete state energy levels. The details are presented here, with respect to 
the ground state. 
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Our bounding procedure is based on the Eigenvalue Moment Method 
(EMM) formahsm of Handy et al (1985,1988a,b). This, Hnear programming 
based (Chvatal (1983)), formahsm has been shown to be exceptionahy weh 
suited for singular perturbation/strong coupling problems. It is very sim- 
ple to use, and involves the application of fundamental theorems arising 
from the classic Moment Problem (Shohat and Tamarkin (1963), Akhiezer 
(1965)), as well as theorems pertaining to the signature structure of bosonic 
(ground state) wavefunctions (Reed and Simon (1978)). 

Specifically, the multidimensional bosonic ground state wavefunction 
must be of uniform signature, which can be taken to be positive: 

^gri-r") > 0. (1) 

It will then satisfy the positive integral relations: 

dxdydz(Vc{'r')y^gr('r') > 0, (2) 

where Vc = ^ „ Ci^rn,nX^y^z^ is an arbitrary polynomial. In terms 
of the power moments, m(pi,P2,P3) = J J J dxdydz x^^y^^z^^ ^^r(^)? 
for nonnegative integer p^'s, these integrals become the Hankel-Hadamard 
(HH), quadratic form, inequalities: 




^ ^ Qi.mi.niMC^i + ^2,mi +m2,ni + n2)Q2,^2,n2 > 0, (3) 

h,mi,ni I2,'m2,n2 



for arbitrary Cs (not all identically zero). 

The Fourier transform of the Schrodinger equation usually admits a 
power series expansion, whose coefficients (i.e. the moments) satisfy a linear 
recursion relation, referred to as the Moment Equation (ME). This relation 
exists for any energy parameter value, E. The entire set of power moments 
is divided into two subsets: 

{fi{t)\ V^} = {M"?)l^ e Ms} U{M^)I^ ^ Ms}. (4) 

The first subset corresponds to the initialization moments, or missing mo- 
ments, which must be specified before all of the other moments can be 
generated, through the ME relation. 

The generated moments, those in the second subset, are linearly de- 
pendent on the missing moments. We can represent the ME relationship 
as 

fiit)= Yl MEit,~t)fiCt), (5) 

~^ ^ M-s- The Me coefficients are dependent on and can be de- 
fined so that the above is also valid for the missing moments as well (i.e. 
ME(ti, S) = ^, for 7i,2 e Ms). 

For one dimensional systems, the number of missing moments is finite, 
and denoted 1 + m^. For multidimensional systems, the number of missing 
moments is infinite; however, at any point in the calculation, one works 
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with a finite number of them, which in turn determine a finite number of 
the generated moments. Generally, there will be many times more generated 
moments, than the corresponding number of missing moments. 

Since the Moment Equation is a homogeneous relation, one must im- 
pose a normalization condition. This is normally done with respect to the 
missing moments. For instance, we can take 

EM7) = 1- (6) 

Assuming this, one can substitute the ME relation, for the generated mo- 
ments, into the HH inequalities. Since all the moments are linear in the 
missing moments, a linear programming problem is defined of the form 

EA-^(£;;C)Mt)>0, (7) 

where the A coefficients are nonlinear ly dependent on E, and quadratically 
dependent on the (arbitrary) Cs. 

In practice, we work within a finite dimension (/) subspace for the 
C coefficients. Define this by {l,m,n) G C^^\ The required moments are 
= ihymi.ni) + (/2,"^2,n2), where (/i,2, ^1,2, ^1,2) G C^^^. One 
must then determine the missing moments that generate these. They in 
turn define the linear programming variable space in Eq.(7). 

For a given dimension, /, at an arbitrary energy value, E, the HH 
inequalities will either have a missing moment solution set, or not 
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U^^^ = (Z). If there is a solution set, it must be convex. This convex set 
may be considered as the intersection of an (uncountably) infinite number 
of polytopes (convex sets formed from the intersection of a finite number of 
hyperplanes) . 

The objective of the linear programming based, algorithmic implemen- 
tation of EMM, is to quickly determine the existence or nonexistence of 
U^E^. At any order /, the feasible energy values (those for which the con- 
vex set exists) define an energy interval, (-E^^'* , E^'' ) , within which the true 
ground state value, Egr, must lie. As the order is increased, the energy 
endpoints define the converging lower and upper bounds to the physical 
answer: 

4" < <...<£,,<... < Eij^» < Eij> , (8) 

/ — ^ OO. 

The EMM formalism was used to generate rapidly converging bounds 
to the ground state binding energy for hydrogenic atoms in superstrong 
magnetic fields (Handy et al (1988a,b)), otherwize known as the Quadratic 
Zeeman effect. This problem had been notoriously difficult, yielding varying 
results depending on the method used. The ability of EMM to define tight 
bounds to the ground state binding energy enabled one to discriminate be- 
tween competing (energy estimation) methods. In particular, it confirmed 
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the results of LeGuillou and Zinn- Justin (1983), which were based on order 
dependent, conformal analysis. 

The consistency of the EMM generated results, and those based on 
conformal analysis, is more than a coincidence. In one dimension, an afSne 
map transformation of the point x is defined by x — ^ , where s and r are 
scale and translation parameters, respectively. An affine map transform of a 
given function, V{x) ^(^7^), corresponds to a translation and stretching 
(or contraction). 

The variational procedure inherent to EMM, is, in fact, affine invariant. 
This is immediately clear from Eq.(2), since the variation samples over all 
polynomial functions, and the space of polynomials is invariant under afSne 
transforms. To this extent, EMM is in keeping with the underlying philos- 
ophy of conformal analysis, and should yield either consistent, or better, 
results. 

This afSne map invariance underscores the fundamental complementar- 
ity between Moment Quantization methods, such as EMM, and explicitly 
multiscale methods such as Wavelet Transform theory (Handy and Murenzi 
(1998)). This further confirms the relevancy of EMM to singular pertur- 
bation type problems which require a careful balancing of large and small 
scale contributions. 

Despite the numerous types of problems the EMM formalism has been 
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applied to, as reviewed in the cited references, it has not been used on prob- 
lems defined on a compact domain. Such problems require a modification 
of the basic formalism, in order to adapt them to the required boundary 
conditions. The infinite quantum lens potential is one such important case. 

Before applying the EMM analysis to this problem, we discuss certain 
preliminary, pedagogic, examples, in order to facilitate the more compre- 
hensive analysis that follows. We provide a short review of EMM by con- 
sidering two, relatively simple, problems. The first is the sextic anharmonic 
oscillator problem, in one dimension. The second is the infinite square well 
potential. These will introduce us to the necessary linear programming, 
HH, relations for Hamburger, Stieltjes, and Hausdorff moment problems. 
The infinite square well problem is an example of the latter, which also 
corresponds to the infinite lens potential. We also discuss how to ellimi- 
nate possible boundary term contributions to the Moment Equation, and 
still be consistent with EMM theory. This analysis has not been presented 
elsewhere, and defines an important, new contribution, with respect to the 
EMM formafism. 
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II. The Sextic Anharmonic Oscillator 



Hamburger Moment Formulation 



Consider the sextic anharmonic oscillator potential problem: 



(9) 



where the kinetic energy perturbation parameter, e, is explicitly noted, for 
later reference. The mass and coupling strength parameters are denoted by 
m and g, respectively. 

The signature structure for the ground (^o) and first excited (^i) 
states are known a priori: ^i(x) = x^Ti(x), where T^(x) > 0. For simplic- 
ity, we confine our analysis to the ground state case. 

Define the Hamburger power moments 



p > 0. Upon multiplying both sides of the Schrodinger equation by x^, 
and performing the necessary integration by parts, we obtain the Moment 
Equation (ME) 




dx x^^(x) 



(10) 



gii{p + 6) = -mii{p + 2) + Eii{p) + e p{p - 1)ijl{p - 2), (11) 

10 



for p > 0. This corresponds to an effective sixth order finite difference 
equation, in which specification of the "initialization" moments, or missing 
moments, {/j,{i)\0 < i < 5}, as well as the energy parameter, E, generates 
all of the remaining moments. 

One important aspect about working within a moments' representa- 
tion is that kinetic energy expansions become regular (i.e. e expansions). 
This is not the case in configuration space, requiring the use of singular 
perturbation type methods (Bender and Orszag (1978)). One immediate 
impact of the regularity in e is that the order of the moment equation does 
not change for e = and e = 0"*" (unlike in configuration space, where the 
order of the differential equation abruptly changes from to 2). 

Let us denote the missing moment order by l + rUs, where rus = 5. We 
can express the linear dependence of the moments on the missing moments 
through the expression 

li{p) = J^MEip,i)li{i), (12) 

i=0 

where 

ME{ii,i2) = Si„i,, (13) 

for < ^1,^2 < ^s- The Me coefficients are readily obtainable, since 
they satisfy the ME relation with respect to the p-index, in addition to the 
preceeding initialization condtions. 
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We must also impose some, convenient, normalization condition. This 
can be chosen to be 

EfW = i- (14) 

Constraining the zeroth order moment, //(O) = 1 — YlT=i 1^^^) -i^^ redefine 
the moment-missing moment relation as 

fiip) = ^MEip,i)m, (15) 

i=0 



where 



and 



M f)^ f Me{p,0), for 1 = 0, , . 

^"1m^(p,£)-M£;(p,0), forl<i<ms- ^ 

From the Moment Problem, we know that the moments of a nonneg- 
ative measure, on the entire real axis, must satisfy the Hankel-Hadamard 
(HH) constraints 

dx (^Cix'\ ^{x) >0, (18) 

for arbitrary C^'s (not all zero), and < / < oo. The zero equality is only 
possible for configurations made up of a finite number of Dirac distributions. 

The HH integral constraints can be transformed into the quadratic 
form expression 
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Ci,fiiii+i2)Ci^ > 0. 



(19) 



ii,i2=0 



These inequalities do not guarantee uniqueness for ^ (i .e. that the physical 
solution is the only one with these moments); however, because we are 
implicitly working with the moments of a physical system, for which there is 
uniqueness, the nature of the Me{p,£) matrix coefficients should guarantee 
uniqueness as well, within the moments' representation, i.e. satisfaction of 
the Carlemann conditions, etc (for further details see Bender and Orszag 



which generates an uncountable number of linear inequalities (i.e. one linear 
inequality for each C-tuple) in the (unconstrained) missing moment variable 
space: 



(1978)). 



One can then substitute the moment-missing moment relation 




I 



(20) 



i=0 ii,i2=0 




(21) 



£=1 



where 




(22) 



n,i2=o 
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and 

I 

B[C]= J2 Ci,MEiii+i2,0)Ci,. (23) 

U,i2=0 

We recall that the missing moments are restricted to ^^^q fJ^ii) = 1- 

Let U^^ denote the (convex) solution set to the above set of HH in- 
equalities, for given E and /. The objective is to determine the feasible 
energy interval, to order /, for which convex solution sets exists: 

£;g(4'\4'))z/4'V0. (24) 

This can be done through a linear programming based cutting method that 
finds the optimal Cs leading to a quick assesment on the existence or non- 
existence ofU^^ (Handy et al (1988a,b)). 

The preceeding formalism is appropriate if the Schrodinger equation 
potential is not symmetric. In the present case, since the potential is sym- 
metric, we can define a more efficient representation by working in terms of 
a Stieltjes moment formulation. This is done in the following section. 

Stieltjes Moment Formulation 

The parity invariant nature of the sextic anharmonic oscillator requires 
that the ground state be symmetric, ^{—x) = ^(x). This in turn introduces 
more moment constraints. 

For symmetric configurations, the odd order Hamburger moments are 
zero, ii{odd) = 0. The even order Hamburger moments can be regarded as 
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the moments of a Stieltjes measure restricted to the nonnegative real axis 
(through a change of variables, y = x^) 



u{p) = M(2p), (25) 

where 

/•OO 

u{p) = / dy yP^y), (26) 
Jo 

and 

The Stieltjes moments also satisfy a moment equation (e = 1): 

gu{p + 3) = -mu{p + 1) + Eu{p) + 2pi2p - 1) u{p - 1), (28) 

P > 0. 

The order of this finite difference moment equation is 1 + = 3, 
leading to the representation 

nig 

u{p) = Y,MEip,i)u{i), (29) 

£=0 

where the u{£) are defined as before, and satisfy the constraint, Ylil^o ^(^) — 
1. 

One important aspect of working with Stieltjes moments is that because 
the underlying function must be positive, all the Stieltjes moments must also 
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be positive (which is not the case for the Hamburger moments). Thus, for 
the adopted normahzation condition, we have 



< u{i) < 1 



(30) 



for < £ < msi= 2). 

Since one is working on the nonnegative real axis, y > , more HH 
constraints are possible. The constraints in Eq.(19), arising from the in- 
tegral expression in Eq.(18), define the necessary and sufficient conditions 
for the moments to correspond to a nonnegative measure on the entire real 
axis. If we pretend that ^{y) exists on the entire real axis, but we want 
it to be zero on the negative real axis, then one must also introduce the 
counterpart to Eq.(18) for the configuration y^{y): 



for cr = 0, 1, and I < oo. Thus, the only way both $(?/) and y^{y) can 
be nonnegative on the entire y-axis is for $(y) = 0, for y < 0. This is 
an intuitive way of motivating the HH-Stieltjes moment conditions for a 
nonnegative measure defined on the nonnegative real axis. Consequently, 
in terms of a quadratic form expression, we have 






(32) 



^1,22=0 
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for cr = 0, 1 and / > 0. 

Repeating the same analysis presented earlier (i.e. substituting the 
moment-missing moment relations, and implementing the linear program- 
ming based cutting procedure) allows us to generate very tight bounds for 
the ground state energy. In particular, for e = m = = 1, one obtains 

1.4356246190092 <E < 1.4356246190178, (33) 

for / = 15. 
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III. Defining Quantizable EMM-Moment Equations 

We now focus on issues of relevance to the application of EMM to 
the quantum lens problem. Consider the configuration F(x) = x^^(x). 
Its Stieltjes moments (for symmetric solutions) will satisfy the moment 
equation derived from Eq.(28) : 

gw{p + 2) = -mw{p) + Ew{p - 1) + 2p(2p - l)w{p - 2), (34) 

p > 2, where w{p) = u{p + 1). This corresponds to an effective 1 + = 4 
order relation since the missing moments {w{0),w{l),w{2),w{S)} must be 
specified before all the other moments can be generated. 

However, application of EMM, to the above moment equation, will not 
generate any discrete state energy bounds. The principal reason for this 
is that the same moment equation ensues if we multiply both sides of (the 
modified Schrodinger equation) 

( - 5^ + mx'^ + gx^ - E^^{x) = V{x), (35) 

by x^'^^, p > 0, provided V{x) is a (symmetric) distribution which is pro- 
jected out when multiplied by x^. Thus, we can have T>{x) = A6{x) + 
B6"{x), where A and B are arbitrary. It is reasonable to expect that 
Eq.(35) admits many bounded, positive, solutions, for arbitrary E; thereby 
explaining the lack of any EMM generated bounds for the w moments. 
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In general, when generating a moment equation, we are free to multiply 
both sides of the the Schrodinger equation by expressions of the form x^T(x) 
(where p > 0) so long as all the zeroes of T are zeroes of the desired physical 
solution (i.e. if T{xz) = 0, then ^(rr^) = 0). If this is not satisfied, then the 
resulting ME relation will fail to distinguish between the true Schrodinger 
equation, and that modified by additional distribution terms supported at 
zeroes of T. 

In accordance with the above, whereas T{x) = x"^ generates a moment 
equation that yields no discrete states, the function T{x) = 1 + does 
generate the ground state solution. Applying x'^p{1 + rr^) to both sides of 
Eq.(9), we obtain the Stieltjes moment equation 

gu{p + 4) = -gu{p + 3) - mu{p + 2) + {E - m)u{p + 1) 

+ [E + 2(p + l)(2p + l)]u{p) + 2p(2p - l)u{p - 1). (36) 

This is a 1 + rus = 4 order relation. Application of EMM generates the 
ground state energy (although at a slower convergence rate): 1.4356178 < 
Egr < 1.4356185, utilizing Stieltjes moments {u{< 30)} . 

A more instructive example is that of the first excited state for the sex- 
tic anharmonic oscillator. The wavefunction will be of the form "^excix) = 
xTi{x), where Ti(x) > 0, and Ti{—x) = Ti(x), for x & We can 
transform the Schrodinger equation into an equation for Ti(x): 

-e(^T;(x) + T'/(x)) + [mx^ + gx^]ri{x) = ETiix). (37) 
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Integrating both sides with respect to will yield the corresponding 
Stieltjes moment equation; however, it will involve (for p = 0) the non- 
(Ti)moment expression dx'^^^^, which is finite. Although a corre- 
sponding EMM analysis can be implemented, it will require a modification 
of the conventional EMM formalism, as previously defined. 

An alternate approach is to simply take T{x) = x, and work with 
the configuration = x'^{x) = rr^Ti(x). The Stieltjes-S moments are 
w{p) = dx x'^^S{x), for p > 0. In terms of the Hamburger moments, 
these become w{p) = fj,{2p-\-l). If we return to the Hamburger ME relation 
in Eq.(ll), and take p = 2p-|- 1, we obtain the desired w-Stieltjes equation: 

-e2p(2p + l)w{p - 1) + mw{p + 1) + gw{p + 3) = Ew{p), (38) 

p > 0. For the case e — m = g = 1, working with the first 30 Stieltjes 
moments, we obtain the bound 

5.033395937697 < Ei < 5.033395937709. (39) 

For the ground state wavefunction, the function T{x) cannot be zero 
except where ^^grix) is zero. For problems defined on a compact domain, 
this means that T can be zero only at the boundary, where the ground 
state wavefunction will, generally, be zero. We discuss this in the following 
section. 
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IV. A Hausdorff Moment Problem: The Infinite Square Well 

We now consider the infinite square well problem 



-dl'^ix) = E^{x), (40) 

where ^(ztL) = 0. The Hamburger moments are ii{p) = J^^ dx xP'^{x). 
For symmetric configurations, we have jii^p) = u{p) = dy yP^{y), 
where y = x'^ and ^{y) = terms of these Hausdorff moments, the 

corresponding moment equation becomes 



-2p{2p - l)u{p - 1) - 2L^P^'{L) = E u{p), (41) 

p > 0. It involves the boundary terms at ±L. 

Relative to the Stieltjes problem, the Hausdorff moment problem intro- 
duces more constraints to the previous Stieltjes (HH) inequalities. We can, 
intuitively, derive these by assuming that ^{y) is nonnegative on [0, (X)). 
This is what the Stieltjes (HH) constraints in Eq.(31) guarantee (i.e. if we 
pretend that the Hausdorff moments are actually Stieltjes moments). 

In order to further constrain such a function so that it be zero on 
the interval [L^, oo), we must require that (L^ — y)^{y) be nonnegative on 
[0,oo): 
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2 

dy{L^-y){J2^iy') ^(^)>0' (42) 

i=0 

I <oo. That is: 

Ci,(L^u{ii + i2)-u{l + ii+i2)^Ci,>0. (43) 

ii,i2=0 

This is the third set of HH constraints that must be added to those in 
Eq.(31), for the Hausdorff problem. 

We can summarize all the Hausdorff-HH relations by 



1 

E + ^2) + r^M^ + h + C\l^ > 0, (44) 

11,12=0 

for cr = 0, 1, 2, where 

r 1, (7 = 

rl') = \0, a=l , (45) 

U^, (7 = 2 

and 

rO, (7 = 

rS^) = J 1, a = 1 , (46) 
l-l, (7 = 2 

for all nontrivial C'^^-^'s, and / > 0. 

We outline how the above constraints lead to the quantization of the 
ground state. 

Let L = 1, and ^ = 2*'(L). Then u(0) = -| > 0, and all the 
remaining moments can be generated, once A is normalized. The first three 
Hausdorff-HH conditions (/ = 0) become: 
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( u{Q) > 0; 

Hausdorff - HH relations {I = 0) ^ < > 0; . (47) 

[u{0)-u{l) >0 

Combining Eu{l) = -2u{0) - A, and u{l) > 0, yields -2^ + u(0) > 0, 
or -| < 1. The third inequahty, 4^ < 1, yields --| + 1<1, orE>0; 



Having established the positivity of E, we are free to impose the normaliza- 
tion A = —1, hence u{0) = ^. Thus, the ME relation effectively becomes 
a zero missing moment problem, with rUs = 0. We can procede with a 
numerical determination of the ground state energy. 

For problems corresponding to = 0, we do not have to implement 
the linear programming based, EMM, formulation. Instead, we can work 
with the nonlinear HH inequalities (which are the relations usually cited in 
the literature) corresponding to the quadratic form relations given previ- 
ously. That is, the Hausdorff-HH linear (in the moments) constraints, are 
equivalent to the nonlinear (in the moments) determinantal relations: 



hence the lower bound 



2<E. 



(48) 




(49) 



where the various HH matrices are 
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KLi. = + ^2) + r'^^^uii + ii + i2), (50) 

for cr = 0, 1, 2, and < ii,i2 < 

The numerical evaluation of these inequalities yields the bounds 

2.4674010541 <E < 2.4674011008, (51) 

utilizing all the HH determinants corresponding to the first seven moments: 
{u{< 6)}. This compares exceptionally well (up to seven decimal places) 
with the true answer, E = (f )^. 

Moment Equations with No Boundary Terms 

In practice, particularly for multidimensional applications, we prefer 
to work with moment equations that do not involve any boundary terms. 
For the infinite square well case, we can do so by multiplying both sides of 
the corresponding Schrodinger equation by xPT{x), where T(±l) = 0. The 
ensuing moment equation will not involve any boundary terms because the 
kinetic energy term becomes 



r+l r+1 

J dxG{x)^"{x) = G'^'\tl-G''^\tl^ J dxG"{x)^{x), 



(52) 



where G{x) = x^T(x). Since both G and ^ are zero at the boundary, no 
boundary terms will contribute to the ensuing ME relation. 
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As an example, let T{x) = 1 — x^, for the L = 1 case. Applying 
— x^) to both sides of the infinite square well problem yields the 
Hausdorff moment equation 

Eu{p +1) = [E- (2p + 2)(2p + l)]u{p) + 2p(2p - l)u{p - 1), (53) 
for p > 0. Application of EMM duplicates the bounds previously cited. 
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V. The Infinite Quantum Lens Problem 

The quantum lens geometry, as shown in Fig. 1, is bounded by the 
z = Ri plane, and the sphere of radius R2- 

Lens Domain = {z > Ri} f^{r < R2}, (54) 

where Ri < R2. 

In a cylindrical coordinate representation, the Schrodinger equation 
for the infinite quantum lens potential problem becomes (in energy units 
Eq = o!t,n2 : and length in units of the radius, a = W Rn — R^ )■ 

- {^-dp{pdp^) + ^dl^ + dl^) = E^ip, 0, z), (55) 

where (r^ = + z^, note that we will be working with the and 
coordinates) 

Rl<r^ < Rl, and Rj < z^ < (56) 
for z > 0. The boundary condition on the wavefunction is 

The radii Ri and R2 can be redefined in terms of the quantum lens 
parameters a and b, where 

a'^ = Rl- Rl and b = R2 - Ri, (58) 
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or, alternatively, 



R2 = and i?i = (59) 



The lens domain transforms into a triangular domain in the {r^,z'^} 
coordinate space, or, equivalently, 

uj = Rl-r'^, and 1^ ^ z'^ - Rl. (60) 

The corresponding domain is 

< a; < and < < - a;. (61) 

Equation (55) is axially symmetric, and the solutions assume the form 
^(p, 0, z) = e-^^<^V^(p, z). 

In the {a;, u} coordinate system, the Schrodinger equation becomes (i.e. 
first transform into {p^,z^} coordinates, then into {r^,z^}, and finally into 



2 

m 



— CO — v 



^{uj,y) = E^{uj,u). (62) 



The boundaries = R2 and 2;^ = R\ become uj = Q and p = Q. 
According to Eq.(61), the {uj^v} physical domain is restricted to the lower 
left triangle of the [0, a^] x [0, a^] square region. The hypotenuse of this 
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triangle corresponds to — uj — u = p"^ = 0. The wavefunction is not zero 
along it; although it is zero along a; = and u = 0. 

Although we shall work within the {uJ,i^} coordinates, in order to derive 
the necessary moment equations, we note that we can rewrite the above 
equation in terms of the coordinates ^ = uj -\- u and rj = uo — v. The 
derivatives become d^j = dn, du = — drj. The p = boundary 
corresponds to ^ = a^. In terms of these new coordinates, the Schrodinger 
equation becomes 

-4([a2 - ^]dl + 2[a2 - ^]d^dr, + [Rj + SUj + - 2r])]d^ - 0^ - 2dr,)^ 



2 

m 



= Ei;. (63) 



We shall refer to the various function coefficients of the derivative operators 
inEq.(63) (i.e. E^,, by 

ra2-^,z = 2,j=0 
2[a2-^],z = l,j = l 

i?i + 3i?? + (^-2r;),i = 0,j = 2. (64) 
-l,i = l,j = 
-2,i = 0,i = l 

Our objective is to derive, for a given quantum number m, a moment 
equation for Eq.(62), involving the moments 

2 2 
na na —UJ 

u{p,q)= / dco dv uj'^i''^'ip{co,iy), (65) 

Jo Jo 

and no boundary terms (in a manner consistent with the requirements de- 
fined in Sec. III). 

28 



In order to achieve the above, for the m = case, we will have to 
multiply both sides of Eq.(62) by G{uj, u) = u^u^, where p,q > 1. We note 
that G{u, ly) = 0, along both boundaries a; = and = 0, where ip = 0. 
Integrating over the triangular domain in {a;, u} does not introduce any 
boundary terms at all, not even along the p = boundary, where ip 0. 
We prove this, below, for each of the contributing terms in the ME relation. 

(i) The terms G{uj, z^)Co,2(4, and G{uj^ ^)Cq^i{^^ r])drj'ip, do not intro- 
duce any boundary terms since those generated by integration by parts (in 
the 7] direction) correspond to points where G{uj, ly) = and ^ = 0. 

(ii) The integration by parts of G{uj,iy)Gi^i{^,r])d^drjip reduces to 

The boundary terms produced by the first term (along the ^ direction) are 
zero since at one point (corresponding to either uj — or u = 0) we have 
G{uj, y) = 0, while at the other (corresponding to ^ = a^), we have Ci,i = 0. 
The boundary terms from the second term are also zero, since at both ends 
(along the rj direction) we have ■0 = 0. 

(in) The integration by parts for G{u, y)[C2,o{C, rj)d'^'ijj — d^tjj] gives us 
d^{GC2,od^ij) - d^{ijd^iGC2,o)) - d^{G^) + [5|(G'C2,o) + d^G]^;. 

The first term introduces no boundary terms (along the ^ direction) because 
at one endpoint we have G{uj, y) = 0, while at the other C2,o = 0. The 
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second and third terms have no boundary term at the point corresponding 
to ip = 0. However, at ^ = a^, since d^C2,o = —1, we obtain a cancellation 
between the only surviving boundary terms. This concludes the proof that 
no boundary terms arise for the m = case. 

For the m case, we take G{uj, v) = uj'Py^{a^ — uj — u), ioi p,q > 1. 
The preceding argument still holds, although the final cancellation of both 
terms is unnecessary because of the additional (a^ — ^) factor introduced 
through the modified G{u^y). 

Given the above, it is now straightforward to generate the required 
moment equation. 
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The m = Moment Equation 



The moment equation for the m = case is 
-f ^(P, q) = Rl P{P - 1) u{p - 2, - [p2 + ^ + 2pq\ u{p -l,q) 

-2Rl pqu{p-l,q-l)^ [q^ + |] u{p, g - 1) + Rlq{q - 1) q-2), (66) 
for g > 1. 

The missing moments {u(0, 0), . . . , u(A^, 0)} and {u(0, 1), . . . , u(0, -/V)}, 
generate all the moments within the square grid [0, A^] x [OjA/"]. We can 
index the missing moments according to Xo = 1^(0,0), Xi = ti(l,0), 
XN = u{N, 0), Xiv+i = uiS), 1), . . ., X2N = 1^(0, iV). 

We can then determine the energy dependent coefficients linking the 
moments to the missing moments 

ms=2N 

u{p,q)= J2 MEip,q,i)Xi- (67) 

The Me coefficients satisfy the moment equation with respect to the p, q 
indices. In addition, ME{pi^,qi^,£2) = Si^^i^, where {pi^,qij denotes the 
coordinates of the missing moments. 

As explained in the previous examples, one can impose a normalization 
condition of the form Xl^=oX^ = 1? constraining xo- Incorporating this 
within the above relation we have 
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ms=2N 

u{p,q) = ME{p,q,0)+ ^ ME{p,q,i)Xi, (68) 

where 

From the positivity theorems of the Moment Problem, we have to im- 
pose the moment constraints arising from the integral relations 

dujdiy Qa{u;,iy) ^ Cijco'^iy^^ ^/j{u;,iy) > 0, (70) 
z,je[o,/]2 

for arbitrary C"s (not all zero), where 

j-l, a = 

K — u — cr = 3 

It is implicitly assumed that ip is zero outside the triangular domain of 
interest. 

These integral inequalities become linear inequalities, with respect to 
the u-moments: 

3 

C'i,,j,(^/^,nW(Ai;^,n+«l +i2,A2;a,n+il +i2))Ci2j2 > 0, (72) 
iiji 12,32 n=l 
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where 



<7,n — 



1,0,0, for a = 
0,1,0, for a = l 
0,0, 1, for (7 = 2 



(73) 



a^, —1, —1 for a = S 



and the A's associated with nonzero fa,nS 




(0,0), for a= 0,n = 1 
(1,0), for a= l,n = 2 
(0,1), for (7= 2,n = 3 



(74) 



(0,0), (1,0), (0,1), forcj = 2> 



If one defines a coordinate pair sequence G [0,/] x [0,/], then 

the set of points covered by (ii^Ji^) + (ihJh) + iM;a,nA2;a,n), he within 
a square grid [0, N]"^, where = 2/ + 1. AU the moments within this grid 
will be generated by the missing moments, previously defined. 

One procedes by substituting the moment - missing moment relation in 
Eq.(68) into Eq.(72). This defines an infinite set of linear inequalities in the 
missing moments, and can be analyzed through the linear programming 
based EMM algorithm. The EMM numerical analysis generates a finite 
number of optimal C"s which determine if, to order /, the (normalized) 
inequalities in Eq.(72) have a solution set, for the specified E value. 
The feasible energies define the converging lower and upper bounds. 

In Tables I and II we give some results of our approach, as a function 
of the ratio ^, for the m = case. Already, for / = 2 (m^ = 10), and 1 = 3 
{rus = 14), we obtain very good bounds for the ground state energy, even 
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for small lens thickness. As the ratio - becomes smaller, the lens becomes 
thiner, with maximum thickness b, and base diameter 2a. 

In Fig. 2, we compare the normalized ground state energy (■^) , for 
m = 0, as a function of the ratio ^, obtained by three methods: 

(z) exact numerical solution (solid line) of Eq.(62); 

(a) perturbation theory (dashed lines), based on a conformal transfor- 
mation into a semi-spherical shape (Rodriguez et al (2001)); 

(iii) EMM analysis, as reported in Table I ( solid black dots for the 
lower and upper bounds, when the "bounding" bars become too small). 

It can be seen that very good agreement is obtained between the exact 
solution and the EMM bounds, for ^ < 0.7. The perturbation results yield 
better agreement with the exact, numerical solution, for ^ > 0.8. 



34 



Table I: Ground state energy bounds (m = 0, / = 2) 



b 


Bounds 


b 

d 


Bounds 


.80 


6 < ^ < 30 


.20 


292.6 < E < 295.5 


.75 


26 < ^ < 33 


.15 


496 < ^ < 506 


.70 


33.8 < E < 35.2 


.10 


1064 < E < 1093 


.65 


38.8 < E <39.4 


.09 


1300 < E < 1340 


.60 


44.33 < E < 44.42 


.08 


1630 < E < 1680 


.55 


50.98< E < 51.03 


.07 


2110 < E < 2180 


.50 


59.54 < E < 59.58 


.06 


2850 < E < 2950 


.45 


70.85 < E < 70.92 


.05 


4070 < E < 4210 


.40 


86.32< E < 86.44 


.04 


6320 < E < 6520 


.35 


108.30< E < 108.70 


.03 


11160 < E < 11520 


.30 


141.60 < E <142.20 


.02 


24950 < E < 25800 


.25 


195.60 < E <196.80 


.01 


99170 < < 102970 
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Table II : Ground state energy bounds (m = 0, / = 3) 



b 
a 


Bounds 


b 
a 


Bounds 


.40 


86.37 <E< 86.39 


.10 


1077 < E <1080 


.30 


141.90 <E < 141.94 


.05 


4120 <E < 4135 


.20 


293.92 <E < 294.00 
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Conclusion 

We have presented a preliminary analysis of the infinite quantum lens 
potential problem in terms of EMM theory. It is anticipated that a different 
choice of coordinates will improve the bounds for | > 0(.7). These con- 
siderations, and the extension of the EMM bounding method to the lowest 
energy state, within each m 7^ class, are the subject of ongoing research, 
and will be communicated in a future work. 
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Caption for Tables and Figures 



Table I: Ground state energy bounds (m = 0,/ = 2). 
Table II: Ground state energy bounds (m = 0, / = 3). 

Fig. 1: Quantum lens geometry of height b and circular coss section of 
radius a. 

Fig. 2: Ground state energy for a quantum lens as a function of the ratio 
^. The energy is given in units of Eq = ^^^2 , as calculated by: a) exact 
numerical solution (solid line); perturbation theory with respect to the ^ 
parameter (dashed line, Rodriguez et al (2001)); and EMM theory, as given 
in Table I (solid black dots). The lower and upper energy bounds are 
represented by "bounding" bars (which cannot be depicted, for smaller ^ 
values) . 
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